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A phase transition of a model compound of the long-range Ising spin glass (SG) 
Dy :c Yi_. c Ru2Si2, where spins interact via the Ruderman-Kittel-Kasuya-Yoshida (RKKY) in- 
teraction, has been investigated. The static and dynamic scaling analyses reveal that the SG 
phase transition in the model magnet belongs to the mean-field universality class. Moreover, 
the characteristic relaxation time in finite magnetic fields, as well as in zero field, exhibits a 
critical divergent behavior, indicating the stability of the SG phase in finite fields. The presence 
of the SG phase transition in the field in this model magnet strongly suggests that the replica 
symmetry is broken in the long-range Ising SG. 
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It has been a great challenge in the statistical physics 
to clarify the nature of glassy states. In this context, 
much attention has been paid to the spin glass (SG) 1 ) 
as the simplest example of glassy systems. According to 
the theoretical studies on the mean-field model of the 
SG, the Sherrington and Kirkpatrick (SK) model, the 
SG state is a thermodynamic equilibrium phase with the 
replica symmetry breaking (RSB). 2 ) The RSB SG state 
is essentially different from uniform states such as param- 
agnetic and (anti)ferromagnetic states where the replica 
symmetry (RS) is preserved. The RSB SG state has a 
complicated multi-valley structure of the free energy, 3 ) 
whereas the RS state has a simple one- or two-valley 
structure. The SK model can explain the glassy nature 
of the SG qualitatively, and the RSB is a vital concept 
for understanding the glassy state. 

However, an application of the SK model to real SG 
systems is not valid as it is, because the interaction range 
is assumed to be infinite in the SK model, whereas it is 
finite in real systems. An alternative theory of the SG, 
the so-called droplet theory based on the short-range in- 
teraction model, 4 ) predicts no RSB and that only two 
thermodynamic states related to each other by a global 
spin flip exist in the SG state, as in uniform states. The 
nature of the SG state in real systems, the RSB state or 
the RS state, is still a controversial issue even at present. 

The most striking difference between the RSB and RS 
SG states is the stability of the SG state under magnetic 
field. In the mean-field picture, a phase transition with 
the RSB occurs on the de Almeida- Thouless (AT) line 
under magnetic field. 5 ) On the other hand, there is no 
SG phase transition at finite temperature under mag- 
netic field in the droplet picture. 4 ) Many numerical and 
experimental studies have been performed to examine 
the stability of the SG state under magnetic field in real 
materials. Detailed studies of an Ising SG with short- 
range interaction Fe^Mni-^TiOa gave evidence against 
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the SG transition in finite magnetic field. 6,7 ) Recent nu- 
merical studies 8 ) support these experimental results and 
claimed that the short-range Ising SG is the RS state 
described by the droplet picture. 9 ) 

However, the nature of the SG state is still an open 
question in other classes, for instance, the long-range in- 
teraction system. The spacial interaction range, as well 
as the spacial dimension, is relevant to the phase tran- 
sition, and hence, the long-range Ising SG could have a 
different nature from the short-range one. The aim of 
this study is to verify whether the thermodynamic SG 
state under magnetic field is preserved in long-range in- 
teraction systems and the possibility of the RSB state 
in real materials. Here, we report on recent experimental 
studies of a long-range RKKY Ising SG Dy x Yi_ 2; Ru2Si2. 
The RKKY interaction via conduction electrons is a rep- 
resentative long-range interaction in real materials. Our 
results clearly indicate the presence of the SG phase in 
magnetic field and strongly suggest the emergence of the 
RSB in the long-range Ising SG, in contrast to the short- 
range Ising SG. 

DyRu2Si2 is a rare-earth intermetallic compound with 
the tetragonal ThCr2Si2-type crystal structure and ex- 
hibits two successive antiferromagnetic transitions with 
T m = 29 K and T N2 = 3.5 K. 10 ) The compound has a 
strong uniaxial magnetic anisotropy and is recognized 
as an Ising antiferromagnet, where Dy 3+ spins inter- 
act via the long-range RKKY interaction. In the di- 
luted compound Dy 2; Yi_ 2 ;Ru2Si2, a SG transition occurs 
with x = 0.103, as shown in Fig. 1(a). Consequently, 
Dy x Yi_ x Ru2Si2 is a good model compound for investi- 
gating the nature of the long-range Ising SG. 

The single crystalline samples of Dy^Yi_ I Ru2Si2 were 
grown by the Czhochralski method with a tetra-arc fur- 
nace. The concentration of Dy was determined by com- 
paring the saturated magnetization of the diluted com- 
pound along the magnetic easy c-axis with that of the 
pure compound DyRu2Si2- Ac and dc magnetizations 
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Fig. 1. (Color online) (a) Phase diagram of Dy a; Yi_ I Ru2Si2- (b) 
T dependences of M/H with the ZFC and FC conditions and 
— X2- (c) T dependences of C mag of the x = 1.0 and 0.103 com- 
pounds. The horizontal and vertical axes are the scaled temper- 
ature and the scaled C ma g, respectively. 

were measured with the SQUID magnetometer MPMS 
(Quantum Design) equipped in the Research Center for 
Low Temperature and Material Sciences in Kyoto Uni- 
versity. The dc magnetization was measured 10 minutes 
after stabilization at a certain temperature and a certain 
field because the samples showed a very slow dynamics, 
as described later. The ac measurements were performed 
with an ac field of 3 Oe and a frequency of 0.01 Hz < 
lu < 1000 Hz. We use a thin plate-shaped sample with 
a weight of 7.0 mg and a size of 5 x 1 x 0.2 mm 3 for 
the dc and ac magnetization measurements to inhibit 
the effects of the diamagnetic field and the eddy cur- 
rent. Specific heat was measured with PPMS (Quantum 
Design). The lattice contribution of the specific heat of 
Dy a; Yi_ a; Ru2Si2 was estimated from the specific heat of 
a non- magnetic reference compound YRu2Si2- 

Figure 1(b) shows the temperature (T) dependence of 
the dc magnetization divided by the magnetic field M/H 
at H = 10 Oe. The dc magnetization was measured under 
the zero-field-cooled (ZFC) and field-cooled (FC) condi- 
tions. A distinct separation of the ZFC and FC mag- 
netizations is found around 1.9 K. In Fig. 1(b), the T 
dependence of the second nonlinear susceptibility xi is 
also shown. X2, which is the order parameter susceptibil- 
ity of the SG, was obtained by fitting the M(H,T)/H 
at each temperature as a function of H 2 in the form of 
M(H, T)/H = xa(T)+ X 2(T)H 2 + • • • . The fitting range 
of H should be sufficiently low to extract X2, and thus, H 
< 50 Oe was chosen to be around 2 K. As shown in Fig. 
1(b), x.2 exhibits a negative divergent behavior toward ~ 
1.9 K. The separation of the ZFC and FC magnetizations 
and the negative divergent behavior of X2 are character- 
istics of the SG phase transition in the zero field limit. 
Figure 1(c) shows the T dependence of the magnetic spe- 
cific heat (C mag ) of the x — 0.103 compound with the 
data of the pure DyRu2Si2- In the plot, the horizontal 
axis is the scaled temperatures T/T g and T/Tm for each 
compound. The SG transition temperature T g of the x = 
0.103 compound is obtained by the static scaling analysis 



described later. The anomaly of C mag of the x = 0.103 
compound at T g is a cusp-like one, being much weaker 
than that of the pure compound at Tni- The cusp- like 
anomaly of C mag at T g is similar to that expected by 
using the SK model; this is in contrast to general SG 
materials, which show no anomaly of C mag at Tg. 1 ) 

To examine the static critical behaviors of the SG 
phase transition in Dy x Yi_ :c Ru2Si2, we analyzed the 
nonlinear susceptibility of the x = 0.103 compound in 
more detail. First, a log-log plot of — X2 against the re- 
duced temperature e = {T/T g — 1) is shown in Fig. 2(a). 
The best plot showing a linear relation between log(— X2) 
and log£ is obtained with T g = 1.925 K and the critical 
exponent 7 = 1.08. As approaching T g closer than e < 
0.1, the log- log dependence of X2 fails. This might be 
the effect of contamination of higher terms. Second, the 
H 2 dependence of the field-dependent nonlinear suscep- 
tibility Xni(if,T) = X o{T) - M(H,T)/H at T g = 1.925 
K is shown in Fig. 2(b) in the double logarithmic scale. 
Xni(H,T) corresponds to the SG order parameter and 
the experimental data exhibits a liner relation in Fig. 
2(b) up to H = 1000 Oe with the critical exponent 5 = 
1.7. Third, we attempted the scaling plot of Xn\(H, T) in 
the form of Xnie - ' 3 vs H 2 e~^ s . As shown in Fig. 2(c), the 
experimental data in the ranges of e < 1.5 and H 2 < 10 6 
Oe 2 collapse on a single curve with (3 = 1.05±0.15 and j36 
= 1.8 ± 0.1, indicating that the nonlinear susceptibility 
exhibits the scaling form of 

Xnl (T,H)=e^F(H 2 E-^). (1) 

The obtained set of critical exponents is roughly sat- 
isfied with the scaling relation of 7 + (3 — (35 and is 
quite different from that in the short-range Ising SG 
Fe^Mni^^TiOs. 11 ) Instead, it is very close to that in the 
SK model. 1 ) 
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Fig. 2. (Color online) (a) Log-log plot of —\2 vs e. (b) Log-log 
plot of Xnl vs H 2 at T = Tg. (c) Full scaling plot of Xnl £_ vs 
H 2 e~ @ s . Dotted lines represent the asymptotic behaviors of the 
scaling function F(x) as x for x — > and x 1 /* for x — > oo. 



Next, we present the dynamical features of the SG 
transition in Dy a; Yi_ ;E Ru2Si2 in zero bias dc field. Fig- 
ures 3(a) and 3(b) show the T dependences of the real 
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Fig. 3. (Color online) (a), (b) T dependences of \' (a) and x" 
(b) in zero bias dc field H = 0. The soild line represents the 
thermal equilibrium susceptibility, (c) Time dependences of q(t) 
at various T values. The horizontal dotted line represents the 
condition for determining t(T). (d) Log-log plot of r vs e. (e) 
and (f) Full scaling plots of qt$l zu vs te zl/ and \"e~^ /xcq vs 
uie~ zu , respectively. Dotted lines in (f) represent the asymptotic 
behaviors of the scaling function G(x) as x for x — > and x$l zv 
for x — » oo. 



and imaginary parts of the ac susceptibility, respectively. 
In Fig. 3(a), the thermodynamic equilibrium susceptibil- 
ity x cq is also shown, which is determined as Xcq(T) = 
dM cq (T,H)/dH \h^o from the dc magnetization under 
the FC condition. The real part of the ac susceptibility x' 
deviates from Xcq at low temperature, and correspond- 
ingly, the imaginary part \" becomes finite. The observed 
dynamics is much slower than that of the conventional 
SGs. 1 ) Firstly, it originates from the long microscopic re- 
laxation time, being estimated at To = 3.7 x 10~ 2 sec 
by the following dynamics scaling analysis. The strong 
Ising anisotropy, and the consequent high energy barrier 
of the spin flipping, is responsible for the long micro- 
scopic relaxation time. A similar slow dynamics was ob- 
served in several Ising magnets. 12, On the other hand, 
we did not find such slow dynamics in a nearly isotropic 
SG Gd a ,Yi_ a; Ru2Si2, which is similar SG material to 
Dy a: Yi_ a ;Ru2Si2 except for its strength of the magnetic 
anisotropy. The enormous slowing down observed near 
T g should be the critical slowing down of the SG phase 
transition, as described later. 

Here, we examine the dynamic scaling analyses of the 
ac susceptibility. According to the dynamic scaling hy- 
pothesis, the characteristic relaxation time r exhibits the 
critical divergence as r = tqE~ Z1 ', where z and v are 
the dynamic and correlation length critical exponents, 



respectively. In this study, we identify t(T) using the 
criterion of x'(l/ r ;7 1 ) = 0.9x C q(T). A similar criterion 
was also used to determine r in Ref. 7. We plot q(t; T) = 
1 — x'(u;T)/xeq(T) with t = l/u>, corresponding to the 
dynamic spin correlation function, against time at sev- 
eral temperatures in Fig. 3(c). The criterion for identify- 
ing t, q(r(T);T) = 0.1, is represented by the horizontal 
dotted line in the figure. Because the frequency depen- 
dence of the ac susceptibility of Dy a; Yi_ 2 ;Ru2Si2 is very 
strong, as shown in Fig. 3(a), we can see the decay of 
the spin correlations over three decades in magnitude at 
4 K in Fig. 3(c). As temperature decreases, the dynam- 
ics becomes slower, and q(t) does not reach 0.1 below 
T g = 1.925 K. This result is consistent with the general 
feature of r, being infinite at T g , and hence, we con- 
cluded that our criterion for identifying r is appropriate 
for the present compound. Figure 3(d) shows a log- log 
plot of t(T) vs e. In this plot, T g was assigned as 1.925 
K, being the same value obtained by the static critical 
scaling analysis, and we obtained the exponent zv = 1.9. 
Furthermore, we show the full dynamic scaling plots of 
q(t;T) and x"( w ;T) in Figs. 3(e) and 3(f), respectively. 
Both q(t; T) and x"( w ! T) in the ranges of e < 2.0 and u 
< 1000 Hz are appropriately scaled with (3 = 0.8 ± 0.1, 
zv = 2.4 ±0.2, and T g = 1.925 K, in the following forms 
of 



q{t;T) = t-^ zv Q{te zv ), 
X>;T)/x eq (T)=^G( We -- ). 



(2) 



The obtained critical exponents are roughly consistent 
with those obtained by the static scaling analysis and 
the analysis of r(T). Surprisingly, the value of zv (w 2) 
is much smaller than that in other SGs, for instance, 
zv « 11 in Fe^Mni^TiOa. 6 ' 7 ) On the other hand, this 
value is consistent with that predicted in the mean-field 
theory, where (zv)mf is 2 with zmf = 4 and vmf = 1/2. ^ 
The static and dynamic critical exponents 7, /3, <5, and 
zv, and the cusp-like anomaly of the specific heat at 
T g indicate that the SG transition in Dy K Yi_ x Ru2Si2 
belongs to the mean-field universality class. Hence, the 
emergence of the RSB is naturally expected. To address 
this, we examine the stability of the SG phase in fi- 
nite magnetic field through the ac susceptibility mea- 
surements. The SG phase transition in magnetic field 
was mostly addressed by observing the irreversibility of 
magnetization; however, it may not be evidence of the 
thermal equilibrium SG phase transition alone because 
of the existence of the nonequilibrium slow dynamics in 
SG materials. 14 - ) A better way to determine the SG tran- 
sition temperature in field T g (H) is the analysis of the T 
dependence of the characteristic relaxation time in field 
r(T, H) to examine whether r(T, H) obeys the dynamic 
critical divergence represented by 



t(T,H)=t (H)\T/T s (H)-1\ 



(3) 



at T > T g (H). If the SG state is stable even in finite 
field, as well as in zero field, the critical divergence of 
r(T, H) is observed. Otherwise, the experimental data 
does not exhibit the critical divergent behavior except 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



' ! 58§8fl 8 . 

x ** S "2-OK 

=0-9Xeq(T) ^Jk 
3.0K 




e=(T/T fl (H)-1) 



0.5 1.5 2.5 3.5 
Temperature (K) 



Fig. 4. (Color online) (a) T dependences of \' in a dc bias field 
H = 500 Oe. The soild line represents the thermal equilibrium 
susceptibility, (b) Time dependences of q(t; T, H) at various T 
values in H = 500 Oe. The horizontal dotted line represents the 
condition for determining r(T, H). (c) Log-log plots of r(T, i?) vs 
£ at various dc bias fields, (d) H dependence of the SG transition 
temperature Tg(H). The inset shows a log-log plot of H vs (1 — 
T g (H)/T g (0)). 



for the case in zero field. In Fe^Mni-^TiOs, the absence 
of SG in finite field was addressed by this way. 6,7 - 1 

Figure 4(a) shows the T dependence of the real part of 
the ac susceptibility x'i^i T, H) at a dc bias field of H = 
500 Oe. The thermal equilibrium susceptibility x eq (T, H) 
= dM cq (T,H)/dH | H=5 ooOc is also shown. In Fig. 4(b), 
q(t = 1/u; T, H) = 1 -%'(„>; T, H)/ X e^(T, H) at H = 500 
Oe is shown. The same criterion for determining t(T, H) 
as in zero field, q(r(T, H); T, H) = 0.1, is adopted, which 
is denoted by the horizontal dotted line in the figure. 
This criterion is employed and the data set of r(T, H) is 
obtained at each dc bias field of H < 1000 Oe. The log- 
log plots of t(T, H) vs T/T s (H) - I at several fields are 
presented in Fig. 4(c). In any field, t(T,H) obeys the 
critical divergent behavior represented by eq. (3) with 
zv w 2.0. The value of zv is almost ii-independent, in- 
dicating that the SG phase transitions in both zero and 
finite fields belong to the same universality class, that is, 
the mean-field one. The H dependence of T g (H) shown 
in Fig. 4(d) exhibits the AT-like behavior, being repre- 
sented as 



H = A[1-T S (H)/T s (0)} a 



(4) 



with a = 1.2 and A = 2500 Oe. The obtained a and 
A are compatible with the calculated values based on 
the SK model of a = 1.5 and A = 3300 Oe. 5 ) These 
results indicate the presence of the SG phase transition 
in finite field and the RSB SG state in Dy a; Yi_ 2 ;Ru2Si2 
in contrast to the short-range Ising SG Fej.Mni_2.TiO3. 



It is not well understood why the SG transition in 
Dy :E Yi_ 2 ;Ru2Si2 belongs to the mean-field universality 

class. The compound belongs to the class of the 3- 

dimensional Ising SG with the long-range RKKY inter- 



action that decays with distance as 



According to 



the theoretical arguments on the d-dimensional SG with 
the long-range interaction represented as r~ a ', the a = d 
systems below the upper critical dimension, correspond- 
ing to the 3-dimensional RKKY interaction systems, be- 
long to the non-mean-field universality class. 15 -* In a sim- 
ilar Ising system, such as LiHo a: Yi_ a ;F4, where the long- 
range dipolar interaction decaying as r~ 3 is dominant, 
the presence of the SG phase is still controversial exper- 
imentally 13, 16 ) and a non-mean-field-type SG transition 
was predicted numerically. 17 ' Further experimental, nu- 
merical and theoretical investigations for the long-range 
SG are required. 

In conclusion, we have investigated the SG phase tran- 
sition in the long-range RKKY Ising SG Dy. E Yi_. E Ru2Si2 
and observed the mean-field critical phenomena. The ex- 
perimental results show the presence of the SG phase 
under finite field and the validity of the SK model in 
Dy_Yi_ x Ru2Si2, strongly suggesting that the RSB SG 
state emerges in real materials with the long-range in- 
teraction as well as in the artificial infinite interaction 
system. 
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